Draft version November 30, 2009 

Preprint typeset using I^Tf^X style emulateapj v. 04/20/08 



On 
O 
O 

> 
O 

o 

CO 

< 

• 

6 : 

(N 
> 
<N 

CN 
<N 

G\ 
O 
On 
O 



DIRECT EVIDENCE FOR OUTFLOW OF METAL-ENRICHED GAS ALONG THE RADIO JETS OF HYDRA A 

C. C. Kirkpatrick 1 , M. Gitti 2 ' 3 ' 4 , K. W. Cavagnolo 1 , B. R. McNamara 1 ' 4 ' 5 , L. P. David 4 , P. E. J. Nulsen 4 , and 

M. W. Wise 6 ' 7 

Draft version November 30, 2009 

ABSTRACT 

Using deep Chandra observations of the Hydra A galaxy cluster, we examine the metallicity structure 
near the central galaxy and along its powerful radio source. We show that the metallicity of the 
intracluster medium is enhanced by up to 0.2 dex along the radio jets and lobes compared to the 
metallicity of the undisturbed gas. The enhancements extend from a radius of 20 kpc from the central 
galaxy to a distance of ~ 120 kpc. We estimate the total iron mass that has been transported out 
of the central galaxy to be between 2xl0 7 M and 7xlO T M which represents 10% - 30% of the 
iron mass within the central galaxy. The energy required to lift this gas is roughly 1% to 5% of the 
total energetic output of the AGN. Evidently, Hydra A's powerful radio source is able to redistribute 
metal-enriched, low entropy gas throughout the core of the galaxy cluster. The short re-enrichment 
timescale < 10 9 yr implies that the metals lost from the central galaxy will be quickly replenished. 

Subject headings: X-rays: galaxies: clusters — galaxies: active — galaxies: abundances 



1. INTRODUCTION 

The elemental abundances of the hot gas at the centers 
of cooling flow cluster s are often enhanced with resp ect to 
those at larger radii (|De Grandi fc Molendil 120011 ). The 
high central iron abundance relative to the a-elements 
suggests the enhancements are largely composed of 
eject a from type la supernova expl osions associated 
with brightest cluster gala xies (BCG) (|De Grandi et all 
120041 : iTamura et all l2004f ). Radial metallicity profiles 
of the hot gas tend to be broader than the stellar 
light profiles of BCGs, indicating that metal-enriched 
gas is diffusi ng outward and mixing w i th the intraclus- 
ter medium ( Rebusco et al. 2005, 2006; iDavid fc Nulsenl 
120081 : iRasera et a l. 2008). BCGs frequently harbor ac- 
tive galactic nuclei (AGN) which are able to heat the 
metal enriched gas surrounding them, causing it to ex- 
pand and mix with low er metallicity gas at larger radii 
([David fc Nul sen 200^). In addition, radio sources asso- 
ciated with AGN may be able to entrain metal enriched 
gas and transport i t outward anisotropically along the 
radio jets and lo bes (|Gopal-Krishna fc Wiitall200ll 
Briiggenl [2Q0l I Omnia et a l. 20041: iHeath et all 
Roediger eraLll2007r >. 

Inclusions of unusually high metallicity gas have 
been found outside the centers of a handful of 
groups and clusters ([Sanders et aTll2004l : IGu et ahl 2007: 
iSimionescu et al1l2008l ). but their association with AGN 
has not been firmly established. Evidence linking a 
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metal enriched outflow to jet activity has been found 
in the Hydra A cluster which shows enhanced metallic- 
ity near its radio lobes i n an XMM-Newton X-ray image 
(Sim ionescu et al"1l2Q09h . Hydra A has an extensive sys- 
tem of cavities embedded in its ICM that is filled with 
radio emission emanati ng from the n ucleus of the BCG 
(iMcNamara et al . 2000; Davi d et all 120011 : iNulsen et all 
[20021 . [20051 : iWise et al.l l2007h . With a total energy of 
over 10 61 erg, Hydra A is experiencing one of the most 
powerful AGN outbursts in the nearby universe. It is 
thus an ideal candidate to search for and characterize a 
metal-enriched outflow. 

We performed such an analysis on a deep, high- 
resolution X-ray image obtained with the Chandra X-ray 
Observatory, and we find clear-cut evidence for metal en- 
riched gas aligned with the radio source. Throughout this 
letter we assume a ACDM cosmology with Ho = 70 km 
s _1 Mpc -1 , Qm = 0-3, and = 0.7. Hydra A is at a 
redshift of z = 0.055, which corresponds to an angular 
scale of 1.07 kpc arcsec -1 . All uncertainties are quoted 
at the 90% confidence level. 

2. DATA ANALYSIS 

2.1. X-ray Data Reduction 

Chandra X-ray observations 4969 and 4970 were pro- 
cessed using CIAO version 4.0.1 and version 3.4.5 of the 
calibration database. Background flares were excluded 
using standard filtering, yielding a net exposure time of 
182 ks. Time-dependent gain and charge transfer ineffi- 
ciency corrections were applied. Blank-sky background 
files used for background subtraction were normalized to 
the source count rate in the 9.5 - 12 keV band. 

2.2. Metallicity Profile Along the Jets 

We have searched for a metal-enriched outflow by com- 
paring the radial distribution of iron along and orthog- 
onal to the radio jets and lobes. Guided by the cavity 
system and 330 MHz radio emission shown in Figure 1, 
we have defined sectors parallel to and roughly orthogo- 
nal to the jets. The northern sector located between 90 
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Fig. 1. — First panel: Residual map of the beta model subtracted surface brightness image from Wise et al. ( 2007). North is towards the 
top, east is towards the left. Second panel: Combined, adaptively smoothed image of Hydra A. The black contour outlines the 330 MHz 
radio emission. The red regions are where the spectra were extracted along the northern jet and the blue regions are where the off-axis 
spectra was extracted. Each bin contains at least 44000 counts. Third panel: Same as the second panel, but for the southern jet. 



and 155 degrees, measured counterclockwise from west, 
was divided into five radial bins with minimum signal-to- 
noise of 210. Two regions east and west of the northern 
jet were created with 45 degree opening angles and bin 
sizes matching those of the northern sector. These re- 
gions are outlined in Figure 1. The spectrum for each 
bin was extracted and corresponding event- wei ghted re- 
spon se matrices were computed. Using XSPEC ([Arnaudl 
1996), single temperature plasma models with absorption 
(WABSxMEKAL) were fit to the spectra in the energy 
range 0.5 to 7 keV. The column densit y was frozen at the 
Gala ctic value of 4.94 x 10 20 cm -2 (jDickev fc Lockmanl 
1990). Temperatures, abundances, and normalizations 
were allowed to vary, while the abunda nce ratios were set 
to the the solar photospheric values of Grevesse & Sauval 
(1998). The spectra are very well fit by this model. 

The projected metallicity profiles along and orthogo- 
nal to the northern jet are shown in the upper panel of 
Figure 2. Between radii of ~ 20 kpc and 120 kpc the 
metal abundance of the hot gas is enhanced by up to 
0.2 dex along the jet relative to the orthogonal sectors. 
Metallicities for the inner and outermost bins are indis- 
tinguishable between the on and off jet sectors. 

The southern sector is located between 230 and 315 
degrees counterclockwise from west. The metallicity pro- 
files along and orthogonal to the southern jet are shown 
in the lower panel of Figure 2. A metallicity enhance- 
ment is also found to the south, extending to approxi- 
mately 120 kpc. The lower points in Figure 2 shown as 
triangles represent the orthogonal sectors and are inter- 
preted as the undisturbed metallicity gradient. The gas 
in the undisturbed regions is warmer by up to 0.5 keV 
compared to the jet regions. 

The metallicity enhancement could be affected by the 
presence of cooler gas. To explore this, we have consid- 
ered a second temperature component when fitting the 
regions along the jets. For example, a component with 
a temperature between 0.5 and 0.7 keV in each sector 
would increase the abundances by approximately 15%, 
but would contribute at most 2% to the total emission. 
However, this and other similar models do not signifi- 
cantly improve the spectral fit. If such an extended cool 



component were present, it would have the effect of in- 
creasing the metal enhancements shown in Figure 2 by 
roughly 0.1 dex. 

2.3. Metallicity Map 

The abundance map shown in Figure 3 was created 
using a weighted Voronoi t essellation binning a lgorithm 
(jCappellari fc Copinl [20031 : iDiehl fc Statlerl 120061 ). Each 
bin has a S/N of 150 and was chosen to balance the 
uncertainties on the best-fit abundances while preserv- 
ing the excellent spatial resolution afforded by the large 
number of source counts. For each bin, a spectrum was 
extracted using the procedure described in Section 2.2. 
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Fig. 2. — Projected metallicity profiles for Hydra A. The circles 
represent regions along the jets and triangles are the average fits to 
the regions east and west of each jet. The dashed lines represent the 
aximuthally-averaged metallicity profile. The left vertical dotted 
lines indicate the outer edges of the cavity pairs A and B from 
Wise et al. (2007), and the right dotted line is for cavity pair C 
and D. 
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Fig. 3. — Metallicity map showing the central 5' x 5' of Hy- 
dra A. Each bin contains approximately 22500 counts. Brighter 
regions represent a higher metallicity. The average error per bin 
is approximately 18%. The 330 MHz emission is shown by the 
white contours and the 1400 MHz emission is shown by the black 
contours. 

The spectrum from each bin was fitted in the energy 
range 0.5-7.0 keV with an absorbed single temperature 
MEKAL thermal plasma model using the parameters de- 
fined previously in Section 2.2. 

The highest metallicity gas, shown in yellow in Figure 
3, extends to the north-east and south of the nucleus. 
This gas is closely aligned with the radio emission and 
is consistent with the metallicity excesses shown in Fig- 
ure 2. The temperature along the north-east jet is ap- 
proximately 2.7 keV, and the southern jet is found to 
be approximately 3.1 keV. The surrounding region has a 
temperature that ranges between 3.4 keV and 3.8 keV. 
T his is similar to what wa s found in the XMM analysis 
bv lSimionescu et'a l. (2009). A plume of uplifted gas with 
a modest metallicity enhancement of ~ 0.1 Z/Z , shown 
in orange, can be seen extending north-west along an X- 
ray filament. Though it is located within what we have 
defined as the undisturbed region described in section 
2.2, it has a negligible effect on our analysis. 

3. DISCUSSION 

3.1. Iron Mass 

We estimate the excess iron mass in two regions. First, 
we assume the metal-enriched gas fills the entirety of the 
on-jet regions defined in Figure 1. Second, we consider 
only in the most significant iron excess shown in yellow 
in Figure 3. We consider both regions because the mass 
estimate depends critically on the volumes and filling fac- 
tors (which we assume to be unity) of the metal-enhanced 
regions. The uncertainty associated with these quanti- 
ties should be smaller in the yellow region alone. At the 
same time, the true enhancement must be larger as it is 
indeed spread over a large volume. We compute the iron 
mass as 

M Fe = pVZf Fe ^ (1) 
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Fig. 4. — The diamond points represent the excess iron mass 
along the northern jet compared to the global profile. The cross 
points represent the excess iron mass along the southern jet. Only 
the second, third, and fourth bins show a significant excess. 



where Z is the metallicity of the gas, f Fe ,o is the iron 
mass fraction of the Sun, V is the volume, and p is the 
density of the gas assuming n e = 1.2nn. The density was 
calculated from the deprojection of the surface brightness 
profile along the jets assuming the volumes of the semi- 
annular wedges shown in Figure 1. 

For the first case the difference between the iron mass 
along the jet and the iron mass of the underlying metal 
distribution in equal volumes is plotted in Figure 4 for 
both north and south regions. A significant excess is seen 
in the second, third, and fourth bins along both jets. The 
total volume of these bins combined is 4.1 x 10 T0 cm 3 . 
We find a total excess iron mass of 3.2^' x 10 7 M 

along the northern jet and 3.6^'g x 10 7 M along the 
southern jet. These values are consiste nt with the upper 
limit found by Si mionescu et al.l (|2QQ9f ). 

In our second case, the volumes of the regions are ap- 
proximated as cylinders extending from 15 kpc outside 
the nucleus to a radius of 100 kpc, with an approximate 
width of 25 kpc. Their total volume is 2.7 x 10 69 cm 3 . Af- 
ter subtracting the average metallicity of the undisturbed 
background 0.42 Z/Z , we estimate the excess iron mass 
within the yellow region in Figure 3 to be ~ 1.6 x 10 7 
M . 

The excesses iron mass represents a significant frac- 
tion of the total iron mass in the gas surrounding the 
BCG. The total iron mass measured within a 30 kpc 
radius is 1.61+0.52 x 10 8 M , consistent with the value 
in both lDavid et all (|2001f ) and lSimionescu et all (|2QQ9h . 
Assuming the the excess iron originated at the center of 
the cluster, this implies that ~ 10% - 30% of the iron 
originating from this central region has been displaced 
along the jets. 

Using the relationship of Bohringe r et al.l (|2004l ). we 
estimate the time required to replenish the lost iron 
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through SNe la and stellar mass loss as, 

t em = (lO- 12 ^ Fe + 2.5 x lO-^TFe)" 1 . (2) 

Here, np e = 0.7 M is the iron yield from SN la and 
7Fe = 2.8 x 10 -3 is the iron mass fraction from stellar 
mass los s. With a ij?-ban d luminosity of Lb = 9.2 x 
10 10 L (David et al. 2001) and a supernova rate of S = 
0.15 SNU (jCappellaro et al] [l999). the enrichment time 
required to replace the uplifted iron is approximately 0.2 
- 0.7 Gyr. The enrichment time is a small fraction of the 
age of the cluster and is only two to seven times larger 
than the age of the outburst that created the surrounding 
shock front. Thus, the central metallicity peak would 
recover if the AGN becomes dormant for a period of only 
a few times its current age. 

Because some star formation is occurring in the BCG, 
we have also considered the contribution from SN II en- 
richment. Assuming one supernova per 100 solar masses 
of sta rs produced and a star formation rate of 1 - 5 M 
yr -1 (jMcNamaralll995D in the BCG, we find only 3% of 
the total iron mass is expected to originate from SN II 
enrichment over this time-scale. 

3.2. Outflow Energy 

The fraction of the AGN outburst energy required 
to lift the enriched gas to its present position pro- 
vides a lower limit to the AGN energy deposited into 
the high metallicity gas near the BCG. This value can 
be estimated by calculating the difference in gravita- 
tional potential energy between the original position of 
the uplifted gas and its current position. Following 
iRevnolds et al.l (|2008[ ), we calculate this quantity as, 



where the sound speed is c s « 750 km s _1 , pi and pf are 
the initial and final densities of the surrounding ICM and 
7 = 5/3 is the ratio of specific heat capacities. Assuming 
all of the displaced gas was initially at the center of the 
cluster, the energy required to lift the metal-enriched gas 
lies between 0.5 - 5xl0 59 erg. This figure is comparable 
to the work required to in flate the inner cav ities against 
the surrounding pressure (|Wise et al.l 120071 ) and is be- 
tween 1% and 5% of the total energy expended by all of 
the cavities and the surrounding shock front. This anal- 
ysis suggests that lo wer power AGN outbursts ty pically 
observed in clusters ([McNamara fc Nulsenll2007t ) would 
have a minor impact on their metallicity distributions. 



4. CONCLUSIONS 

Hydra A has undergone a major AGN outburst that 
has transported 2 - 7xl0 7 M of metal-enriched gas out 
of the BCG along the radio jets. The gas has been en- 
hanced by up to 0.2 dex in metallicity compared to the 
ambient gas out to distances as large as 120 kpc. Ap- 
proximately 10 59 erg was expended while uplifting the 
metal-enriched gas, which is 1% to 5% of the total AGN 
outburst energy. 



Met 



In 



(3) 



This work was supported by generous grants from the 
Natural Sciences and Engineering Research Council of 
Canada and Chandra grant G07-8122X. M. Gitti ac- 
knowledges support by grant ASI-INAF 1/088/06/0. 



REFERENCES 



Arnaud, K. A. 1996, Astronomical Data Analysis Software and 

Systems V, 101, 17 
Bohringer, H., Matsushita, K., Churazov, E., Finoguenov, A., Sz 

Ikebe, Y. 2004, A&A, 416, L21 
Briiggen, M. 2002, ApJ, 571, L13 
Cappellari, M., Sz Copin, Y. 2003, MNRAS, 342, 345 
Cappellaro, E., Evans, R., Sl Turatto, M. 1999, A&A, 351, 459 
David, L. P., Nulsen, P. E. J., McNamara, B. R., Forman, W., 

Jones, C, Ponman, T., Robertson, B., & Wise, M. 2001, ApJ, 

557, 546 

David, L. P., Sz Nulsen, P. E. J. 2008, ApJ, 689, 837 
De Grandi, S., & Molendi, S. 2001, ApJ, 551, 153 
De Grandi, S., Ettori, S., Longhetti, M., & Molendi, S. 2004, 
A&A, 419, 7 

Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 
Diehl, S., Sz Statler, T. S. 2006, MNRAS, 368, 497 
Grevesse, N., Sz Sauval, A. J. 1998, Space Science Reviews, 85, 
161 

Gopal-Krishna, Sz Wiita, P. J. 2001, ApJ, 560, L115 
Gopal-Krishna, Sz Wiita, P. J. 2003, Radio Astronomy at the 
Fringe, 300, 293 

Gu, J., Xu, H., Gu, L., An, T., Wang, Y., Zhang, Z., & Wu, X.-P. 

2007, ApJ, 659, 275 
Heath, D., Krause, M., & Alexander, P. 2007, MNRAS, 374, 787 
McNamara, B. R. 1995, ApJ, 443, 77 
McNamara, B. R., et al. 2000, ApJ, 534, L135 
McNamara, B. R., Sz Nulsen, P. E. J. 2007, ARA&A, 45, 117 
Nulsen, P. E. J., David, L. P., McNamara, B. R., Jones, C, 

Forman, W. R., & Wise, M. 2002, ApJ, 568, 163 



Nulsen, P. E. J., McNamara, B. R., Wise, M. W., & David, L. P. 

2005, ApJ, 628, 629 
Omma, H., Binney, J., Bryan, G., Sz Slyz, A. 2004, MNRAS, 348, 

1105 

Rasera, Y., Lynch, B., Srivastava, K., Sz Chandran, B. 2008, ApJ, 
689, 825 

Rebusco, P., Churazov, E., Bohringer, H., Sz Forman, W. 2005, 

MNRAS, 359, 1041 
Rebusco, P., Churazov, E., Bohringer, H., Sz Forman, W. 2006, 

MNRAS, 372, 1840 
Reynolds, C. S., Casper, E. A., & Heinz, S. 2008, ApJ, 679, 1181 
Roediger, E., Briiggen, M., Rebusco, P., Bohringer, H., Sz 

Churazov, E. 2007, MNRAS, 375, 15 
Sanders, J. S., Fabian, A. C, Allen, S. W., & Schmidt, R. W. 

2004, MNRAS, 349, 952 
Simionescu, A., Werner, N., Finoguenov, A., Bohringer, H., Sz 

Briiggen, M. 2008, A&A, 482, 97 
Simionescu, A., Werner, N., Bohringer, H., Kaastra, J. S., 

Finoguenov, A., Briiggen, M., Sz Nulsen, P. E. J. 2009, ASzA, 

493, 409 

Tamura, T., Kaastra, J. S., den Herder, J. W. A., Bleeker, 
J. A. M., Sz Peterson, J. R. 2004, A&A, 420, 135 

Wise, M. W., McNamara, B. R., Nulsen, P. E. J., Houck, J. C, Sz 
David, L. P. 2007, ApJ, 659, 1153 



